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Abstract 

Most existing work on adaptive allocation of subcarriers and power in multiuser orthogonal 
frequency division multiplexing (OFDM) systems has focused on homogeneous traffic consisting 
solely of either delay-constrained data (guaranteed service) or non-delay-constrained data (best-effort 
service). In this paper, we investigate the resource allocation problem in a heterogeneous multiuser 
OFDM system with both delay-constrained (DC) and non-delay-constrained (NDC) traffic. The 
objective is to maximize the sum-rate of all the users with NDC traffic while maintaining guaranteed 
rates for the users with DC traffic under a total transmit power constraint. Through our analysis we 
show that the optimal power allocation over subcarriers follows a multi-level water-filling principle; 
moreover, the valid candidates competing for each subcarrier include only one NDC user but all 
DC users. By converting this combinatorial problem with exponential complexity into a convex 
problem or showing that it can be solved in the dual domain, efficient iterative algorithms are 
proposed to find the optimal solutions. To further reduce the computational cost, a low-complexity 
suboptimal algorithm is also developed. Numerical studies are conducted to evaluate the performance 
the proposed algorithms in terms of service outage probability, achievable transmission rate pairs 
for DC and NDC traffic, and multiuser diversity. 

Index Terms 

Orthogonal frequency division multiplexing (OFDM), constant-rate transmission, variable-rate 
transmission, power control, convex optimization, water filling. 

I. INTRODUCTION 

Future broadband wireless networks are expected to support a wide variety of communica- 
tion services with diverse quality-of- service (QoS) requirements. Applications such as voice 
transmission and real-time video streaming are very delay-sensitive and need guaranteed 
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throughput. On the other hand, applications like file transfer and email services are relatively 
delay tolerant so variable-rate transmission is acceptable. From the physical layer point of 
view, transmission of delay-tolerant or non-delay-constrained (NDC) traffic can be viewed as 
an ergodic capacity problem [1], where maximizing the long-term average transmission rate 
is the goal. Thus, wireless resources, such as transmission power and frequency bandwidth, 
can be dynamically allocated so as to exploit the time or frequency selectivities of broadband 
wireless fading channels. Likewise, transmission of delay-sensitive or delay-constrained (DC) 
traffic can be regarded as a delay-limited capacity problem [2] in which a constant trans- 
mission rate should be maintained with probability one regardless of channel variations. In 
this case, it is desirable to allocate more transmission power and frequency bandwidth when 
the channel experiences deep fade and to allocate less resources when the channel is under 
favorable conditions. We investigate in this work resource allocation in a broadband wireless 
network that supports simultaneous transmission of users with delay differentiated traffic. Our 
focus is on the formulation of an analytical framework from the physical layer perspective 
as well as the design of efficient and practical algorithms. 

Multicarrier transmission in the shape of orthogonal frequency division multiplexing (OFDM) 
is a leading technique to provide spectrally efficient modulation as well as user multiplexing 
in future wireless systems. With OFDM technique, the broadband wireless channel is divided 
into a set of orthogonal narrowband subcarriers. In a single user system, since the channel 
frequency responses are different at different subcarriers, the system performance can be 
significantly enhanced by adapting the transmission parameters such as modulation, coding, 
and power over each subcarrier. For instance, the transmitter can send at higher transmis- 
sion rates over subcarriers with better channel condition while lower rates or no data over 
subcarriers in deep fade. This follows the well-known water-filling principle. In a multiuser 
system, different subcarriers can be allocated to different users to provide a multiple access 
method, also known as OFDMA. As the channels on each subcarrier are likely independent 
for different users, the subcarriers experiencing deep fade for one user may not be in deep 
fade for other users. As a result, each subcarrier could be in good condition for some users 
in a multiuser OFDM system. By adaptively allocating the subcarriers among multiple users 
based on instantaneous channel information, multiuser diversity can be utilized to boost the 
overall system efficiency. 

Adaptive resource allocation in multiuser OFDM systems has focused on homogeneous 
traffic only. In such systems, the traffic consists solely of either DC data requiring constant- 



rate transmission [3], or variable-rate NDC data which can be served in a best-effort man- 
ner [4]-[6]. For systems with pure DC traffic, the problem is to minimize the total transmit 
power while satisfying a basic transmission rate for each user, which is often referred to as 
margin adaption [7]. In [3], an iterative algorithm was proposed to allocate each user a set of 
subcarriers and then determine the power and rate for each user on its allocated subcarriers. 
For systems with pure NDC traffic, the problem is often formulated as maximizing the sum- 
rate of the system subject to a total transmit power constraint. This formulation is also known 
as rate adaptation [7]. In [4], it was shown that the total sum-rate of a multiuser OFDM 
system is maximized when each subcarrier is allocated to the user with the best channel gain 
for that subcarrier. The total transmit power is then distributed over the subcarriers using the 
water-filling algorithm. This result holds, however, only for single-antenna systems. It is no 
longer optimal when multiple antennas are deployed at the base station due to the spatial 
multiplexing gain [5], [6], [8]. Other problem formulations for systems with pure NDC 
traffic take user fairness into account. For example, [9] studied the max-min criterion which 
aims to maximize the transmission rate of the bottleneck user. In [10], it was proposed to 
maintain proportional rates among users for each channel realization. A utility-function based 
optimization framework to balance system efficiency and user fairness was also discussed in 
[11]. 

In this paper, we consider the subcarrier and power allocation problem in a heterogeneous 
multiuser OFDM system where DC and NDC traffic is supported simultaneously. Users in the 
system are classified into DC users and NDC users based on their traffic delay requirements. 
We assume that the total transmit power from the base station is fixed. Our objective is to 
maximize the sum-rate of all the NDC users while maintaining the basic transmission rates 
of all the DC users over every transmission frame. A similar problem was studied in [12]. 
However, it assumed static subcarrier allocation so only the transmit power adaptation was 
discussed. Our work, instead, considers joint subcarrier and power adaptation in multiuser 
OFDM systems, which is one step forward of the previous work. This multiuser subcarrier 
and power allocation problem is a mixed integer programming problem, the complexity of 
which increases exponentially with the number of subcarriers. To make the problem more 
tractable, we transform it into a convex programming problem by introducing time- sharing 
variables. We show that, for a given subcarrier assignment, the optimal power distribution is 
achieved by multi-level water-filling. In particular, the water level of each DC user depends 
explicitly on the channel gains of its assigned subcarriers and its basic rate requirement, 



and can differ from one another. On the other hand, the water levels of all NDC users are 
the same. We also show that, for the optimal subcarrier assignment, the set of valid user 
candidates competing for each subcarrier consists of all the users from the DC group and 
one from the NDC group with the best channel gain. Using these properties, we propose an 
efficient iterative algorithm to compute the optimal solution numerically. Alternatively, the 
original problem is solved in the dual domain by using dual decomposition. It is shown that 
the dual updates can be done efficiently using an ellipsoid algorithm. In addition, we present 
a suboptimal algorithm with linear complexity in the number of subcarriers and the number 
of users. 

The rest of this paper is organized as follows. In Section II we introduce the system model 
and describe the problem formulation. In Section III we formulate the resource allocation 
problem as a convex optimization problem by using time-sharing technique and present 
analytical frameworks of the optimal solution. An iterative algorithm to search the optimal 
solution is also presented. In Section IV, we attempt to solve the problem using dual approach. 
A low-complexity suboptimal algorithm is given in Section V. In Section VI, we present 
numerical results of our proposed algorithms in a multiuser OFDM system. Finally conclusion 
and discussions are given in Section VII. 

II. System Model and Problem Formulation 

We consider the downlink of a multiuser OFDM system with block diagram shown in 
Fig. [TJ The system consists of K mobile users. The first Ki users have DC traffic, which 
requires a constant transmission rate of Rk (k = l,...,Ki) bits per OFDM symbol, re- 
spectively. The traffic of the remaining K — K\ users has no delay constraint and can be 
delivered in a best-effort manner. Note that K is the number of users that are scheduled for 
transmission during a certain transmission interval. The total number of users in a practical 
system may be much larger than K and, hence, other multiple access techniques such as 
time-division multiple access (TDMA) are needed in conjunction with OFDMA. The data 
streams from the K users are serially fed into the encoder block at the base station transmitter. 
The total channel bandwidth is B Hz and is divided into N orthogonal subcarriers, which 
are shared among the K users. The transmission is on a time-frame basis, where each frame 
consists of multiple OFDM symbols. The fading coefficients of all users are assumed to 
remain unchanged within each transmission frame but can vary from one frame to another. 
All channel information is assumed perfectly known at the central controller, which can 



be embedded with the base station. Typically, the channel information can be collected by 
estimating it at each user terminal and sending it to the base station via a feedback channel, 
or through channel estimation of the uplink in a time-division duplex system. Based on the 
instantaneous channel inputs, the central controller allocates different subcarriers to different 
users and determines the amount of power/bits to be transmitted on each subcarrier through 
the subcarrier and power/bit allocation algorithm. The resulting allocation information is used 
to configure the encoder block at the base station transmitter and to facilitate the subcarrier 
selector and decoder at each user receiver. Note that this allocation information may be sent 
to each user via a separate channel. The output data symbols from the encoder are then 
modulated by the inverse fast Fourier transform (IFFT). Guard interval is inserted to ensure 
orthogonality between OFDM symbols. The total transmit power from the base station is 
fixed and is given by P T . 

The broadband wireless channel between the base station and each user terminal is assumed 
to be frequency- selective Rayleigh fading. However, the channel in each subcarrier is narrow 
enough to experience flat fading. Let rk, n denote the transmission rate of user k on subcarrier 
n in bits per OFDM symbol. It depends on the channel gain hk, n an d the allocated power 
P kt n of user k on subcarrier n. In general, r fc n can be expressed as 



where iV is the power spectral density of additive white Gaussian noise and T is a constant, 
usually called the signal-to-noise ratio (SNR) gap [7]. When instantaneous mutual information 
is used to characterize the achievable transmission rate, we have T = 1 (0 dB). If practical 
signal constellations are used, T is a constant related to a given bit-error-rate (BER) require- 
ment. For example, when uncoded QAM modulation is used we have T = — ln(5-BER)/1.5. 
In general, the gap serves as a convenient mechanism for analyzing the difference between 
the SNR needed to achieve a certain data rate for a practical system and the theoretical limit. 
Throughout this paper we use (0Q) as a unified form to characterize both the theoretical mutual 
information and practical transmission rate. 

The problem we consider here is to optimize the allocation of subcarriers and power 
under the total transmit power constraint so as to maximize the sum-rate of all the K — K\ 
NDC users while satisfying the individual rate requirement for each of the K\ DC users. 
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Mathematically, the given problem can be formulated as 
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where Vt k is the set of subcarriers assigned to user k. f2 fc 's must be mutually exclusive 
since each subcarrier is allowed to be used by one user only. In general, it is necessary to 
share the same subcarrier among multiple users in order to achieve the multiuser capacity 
region [13]. This suggests that superposition coding together with high complexity decoding 
should be used. However, there is only a small range of frequency with overlapping sharing 
according to [13] when optimal power control is used. We therefore focus on mutually 
exclusive subcarrier assignment schemes, which can also simplify transmitter and receiver 
implementation for practical systems. 

III. Time-Sharing Based Optimal Subcarrier and Power Allocation 

Finding the optimization variables Vt k an d Pk, n for all k and n in © is a mixed integer 
programming problem. In the system with K users and N subcarriers, there are K N possible 
subcarrier assignments since each subcarrier can be used by one user only. For each subcarrier 
assignment, the total power will be allocated to meet the individual rate requirement for each 
DC user and at the same time to maximize the sum-rate of the NDC users. The subcarrier 
assignment together with its associated power allocation that results in the largest sum-rate 
while satisfying all the constraints is the optimal solution. 

An approach to make the problem more tractable is to relax the constraint that each 
subcarrier is used by one user only. We introduce a sharing factor p k , n £ [0, 1] indicating 
the portion of time that subcarrier n is assigned to user k during each transmission frame. 
This time-sharing technique was first proposed in [3] and has been frequently used in the 
context of subcarrier assignment in multiuser OFDM systems to convert a mixed integer 



programming problem into a convex optimization problem [5], [6], [9], [10]. In addition, 
we introduce a variable Sk >n and define it as s fc n = Pk, n Pk,n for all k and n. Clearly, Sk >n 
becomes the actual amount of power allocated to user k on subcarrier n, whereas Pk >n is the 
power as if subcarrier n is occupied by user k only. If p kn = 0, we always have n = but 
Pfc n is not necessarily equal to zero. For notation brevity, we let a^n = \hk, n \ 2 / (r kN^B / N) 
for all k and n, which is called the effective channel-to-noise ratio (CNR) of user k on 
subcarrier n. Here, for the purpose of generality, the subindex k is added to the SNR gap T 
to include the case when each user has different BER requirements if adaptive modulation 
and coding is used. With the aid of time-sharing factors pk >n we now readily transform the 
original problem © into: 
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The objective function © is a sum of functions of the form f(pk, n , Sk : n) = Pfc,nlog 2 (l + 
Csk, n / Pk,n) i where C is some positive constant. By evaluating the Hessian matrix of f(pk,n, Sk,n) 
at pk, n and Sk, n , we can prove that f(pk,n, Sfc,n) is concave [14]. Thus, the objective function is 
also concave since any positive linear combination of concave functions is concave. Moreover, 
since the inequality constraint functions in © are convex and the constraints in ©-© are 
all affine, the feasible set of this optimization problem is convex. Therefore, the problem in 
©-© is a convex optimization problem and there exists a unique optimal solution, which 
can be obtained in polynomial time. In the following we derive some desirable properties of 
the optimal solution. 



The Lagrangian of the above problem is given by 
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where (3= (0i, . . . , ^ 0, /x > and « = (t>i, . . . , v N ) are the Lagrange multipliers for 
the constraints ©, © and ©, respectively. The boundary constraints © will be absorbed 
in the Karush-Kuhn-Tucker (KKT) conditions [14] as shown later. Let p* kn and s* kn denote 
the optimal solution, if it exists, for 1 < k < K, 1 < n < N. Applying the KKT conditions, 
we can obtain the necessary and sufficient conditions for p* k n and s* k n . Specifically, p* k n and 
s* kn should satisfy the following equations: 
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A. Optimal Power Allocation for Given Subcarrier Assignment 

In this subsection, we present the optimal power distribution when subcarrier assignment 
is given. 

Let {pk, n } be any given subcarrier assignment scheme. Differentiating the Lagrangian in 
© with respect to Sk, n and substituting the result into the KKT condition ©, we obtain: 
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Equations (fT2l) and (TOT) clearly show that the optimal power allocation follows the standard 
water-filling approach, except that the allocated power is only on for p fc n portion of time. 
For each user, more power will be allocated to the subcarriers with higher CNRs and vice 
versa. But the water levels of different users can be different. Specifically, the water level 
of each DC user is given by L k = /3fc/(/iln2), for k = 1, ... , Ki, and it should ensure the 
basic rate requirement R k in ((4]). Substituting dT2]) into the KKT condition (fTT|) and in view 
of Pk 7^ 0, we obtain the closed-form expression for L k given by: 

Lk = [rw Kn) Pfc " 

where Q k is the set of subcarriers that is assigned to user k with p k)U > and satisfies 
«fc,n > I n me case wnere the given subcarrier assignment is mutually exclusive, i.e., 

all the pfc jn 's only take one or zero, the water levels L k can be re-expressed as 

/ o Rk \ 1/gk 
Lk = TT^H ' (15) 

Vll„'=l a k,n'J 

where a^v > ot k ^ >•••,> &k,\n k \> are the ordered CNRs of the A;-th DC user on its allocated 
subcarrier set Q k , and g k (< \£l k \) is the largest integer satisfying a k , g ' > 1/L k . The water 
level of all the NDC users, on the other hand, is observed from (fT3l) to be the same and is 
given by L = l/(pln2). Let Pdc,t = J2k=i E!!=i Pk,nP kn represent the actual total power 
consumed by the K± DC users. Then the water level of NDC users, L , can be obtained 
numerically by using the total power constraint P D c,t + Ef=K 1+ i E«=i Pk,nP k ^ n = Pt- 

Fig. |2] illustrates the optimal power distribution based on multi-level water-filling in a 
multiuser OFDM system with K = 4 users and N = 8 subcarriers. In the figure, the height 
of each blank region represents the inverse of the channel-to-noise ratio and the height of 
shadowed ones is the allocated power. The water levels of the two DC users, k — 1 and 
k = 2, are given by L\ and L 2 , respectively, and they are determined explicitly by the basic 
rate targets. The water level of the two NDC users, k = 3 and k = 4, is given by L , 
and it depends on the total available power after the subtraction of the power consumed by 
DC users. This interpretation on the determination of water levels implicitly imposes higher 
priorities on DC users, for which the basic rate targets must be guaranteed all the time. 

B. Optimal Subcarrier Assignment 

The subcarriers and power should be allocated jointly to achieve the optimal solution of the 
problem formulated in ©-CO). The previous subsection discussed the analytical expressions of 
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the optimal power allocation for a given subcarrier assignment. In this subsection, we derive 
the optimal strategy for subcarrier assignment assuming the power allocation is optimized. 
Taking the partial derivative of the Lagrangian in ® with respect to p k)Tl , we have: 
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where (3 k = fi k for 1 < k < K 1 and j3 k = 1 otherwise. Now we substitute the optimal power 
allocation (fT2)) and (fT3l) into (fT6l) and apply the KKT condition (flOl) . then we get: 

1, v n < Hk, n (L , L k ) 
0, v n > H k , n (Lo, L k ) 

Here, the variable L k , for k — 1, . . . , K is defined as L fc = L fc for 1 < k < K 1 and L k = L 
otherwise, and the function H k:Tl (L , L k ) is given by: 

1 

C-k,nLk t 

Due to the constraint in ©, we conclude that the optimal subcarrier assignment is decomposed 
into N independent problems. That is, for each subcarrier n, if H k:n (L , L fc )'s, for k = 
1,...,K, are all distinct, then only the user with the largest H ktn (L , L k ) can use that 
subcarrier. In other words, we have: 

P*k'(n),n = 1, Pt,n = 0, Vfc ^ k'(n) (18) 

where k'(n) = argmaxi<fc<x H k:n (L , L k ). Hence, it follows that for a given set of water 
levels {L , L 1 , . . . , L Kl }, or equivalently, a given set of Lagrange multipliers {/3,fi}, we 
can determine the optimal subcarrier allocation using (TT8l) . However, the optimal solution 
obtained may not satisfy the individual rate constraint @> and the total power constraint Q. 

Function H k:n (L , L k ) in (flTT) plays a key role in finding the optimal solution of subcarrier 
and power allocation. We now take a closer look at its features. Firstly, by differentiating 
H k! n(L , L k ) with respect to a kjn at a k:1l > 1/L k , it is seen that H ktn (L , L k ) is monotonically 
increasing in a k >n . As a result, for each subcarrier, the user with a larger CNR is more likely 
to be assigned this subcarrier. In the extreme case where DC users are absent (i.e. K\ = 0), 
each subcarrier will be assigned to the user with the largest CNR. This agrees with the 
previous result in [4]. In the general case, the candidate pool for accessing each subcarrier 
consists of all the K\ DC users and the only one NDC user with the largest CNR. This 
consequently suggests that the exhaustive search for optimal subcarrier assignment in the 
original problem © has a complexity of [K\ + 1)^ rather than K N . Secondly, we also 



observe that H kn (L , L k ) is a non-decreasing function of L k and L , respectively, when 
1 < k < K\ and K\ < k < K. This is obtained by differentiating H kn (L , L k ) with respect 
to L k and L for 1 < k < K\ and K\ < k < K, respectively. Thus, increasing the water 
level L k of DC user k increases its chance to occupy more subcarriers. Similarly, increasing 
the water level L of NDC users allow them to access more subcarriers. 

Using the above observations, we propose in the next subsection an iterative algorithm to 
find the optimal water levels and determine the corresponding subcarrier assignment so that 
all the rate and power constraints are satisfied. 

C. Iterative Algorithm 

The thrust of the algorithm is to obtain the set of optimal water levels {L , L±, . . . , Lk x } 
using two nested loops. The outer loop varies L to meet the total power constraint. The 
inner loop searches {L x , . . . , L Kl } and determines the optimal p kjU for all k and n at a given 
value of L to satisfy the basic rate requirement for every DC user. The algorithm is outlined 
as follows. 

Optimal Subcarrier Assignment Algorithm 
Main function 

a) Set L LB = L UB = minii- 1<fc </< {l/a k)Tl } - A 

l<n<N 

b) Find the optimal {L k , p k)U } at L = L UB 
Compute the actual power consumption P' T using (Q3 and (TT3T) 
if Pip < P T , update L UB = 2L UB , and repeat Step b) 
else, go to Step c) 

c) Find the optimal {L k , p Kn } at L = (L LB + L UB )/2 
Compute P' T using CL2]) and C3 
if P' T > P T , let L UB = L 
elseif P' T < P T , let L LB = L 
Repeat Step c) until P' T = P T 

Function: find the optimal {L k , p k>n } at a given L 
1) Set L k using (O with Vt k = {1, . . . , N}, for k = 1, . . . , K x 
Compute H kjn using (flTT) . Vfc, n 
Obtain p kjTl using (TT8T) . Vfc, n 



'A is a very small number. 



\1 

2) Compute R' k = J2n=i pk^og 2 (L k a k ^)} + for k = 1, . . . , K x 
3a) Find k* with R' k , < R k * and R' k , - R k * < R' k - R k for all 1 < k < K x 
3b) Find the subcarrier set A k * = {n\p k *^ n < 1} for the found k* 
for each n £ A k * 

Let = argmaxx^fc^x-fffc.n and obtain L k *(n) such that i?fc* in = (1 + A)H k ( n ) :n 
3c) while i?^.* < R k * for the found fc* 

Denote m = arg min„ e- 4 fc , L k * (n) 
Update L k * = L k * (m) 

Update p fc%m = 1, p Km = for \/k^k*, and A* = A-* - {m} 
Compute Bf k , = Yln=iPk*,n\\og 2 (L k *a k * tn )} + 
3d) if # fc . > i2fc., let i^ fc . = i?^ - log 2 (L fc *a fe . im ) 
if i?^,* < R k * 

Update p fc%m = (i? fc . - R' k *)/\og 2 (L k *a k * )m ) and p fc(m ), m = 1 - p k *,m 
elseif > i? fc . 

Update p k * >m = 0, Pk( m ),m = 1 and compute using (IT4]) 
4) Repeat Steps 2)-3) until i?' fc = i? fc for all k = 1, . . . , K x 

In the outer loop (main function), we initialize L to a value slightly below the minimum 
reciprocal of the CNRs of all NDC users over all subcarriers so that none of the NDC users 
is assigned any power resource based on (fT3l) . By doing this, the subcarrier and power will 
be initially allocated to all DC users as if NDC users were absent. We then keep increasing 
L until the actual total power consumption P' T exceeds the total available power P T and 
an upper bound of L is obtained. Note that if the number of subcarriers in a practical 
system is large enough, we can always find at least one subcarrier fully occupied by one 
NDC user and, therefore, an analytical upper bound of L can be derived from (IT3T) as 
£ub = Pt + niax^ 1< fc<ii'{l/a/ c>n }. The algorithm then proceeds to use the bisection method 
to update L so that P' T converges to Pt- The outer loop converges because the actual total 
power consumption P' T increases monotonically with L given the target data rates for all 
DC users are satisfied. 

The operation in the inner loop (function: find the optimal {L k , p k n } at a given L ) is 
similar to the algorithm introduced in [3]. Each L k is first initialized to the minimum water 
level needed by DC user k to achieve its target data rate, which happens when all the N 
subcarriers in the system are assigned to this user. We then gradually increase L k for one of 



the DC users until the target data rate for this user is satisfied. Increasing L k is carried out 
step by step and the increment at each step is the minimum required value allowing only 
one more subcarrier to be added to this DC user. During this process, there are chances that 
the newly added subcarrier is not fully needed and should be time-shared with other users. It 
is also possible that this subcarrier is not needed at all and it only needs to adjust L k using 
(TT4)) . The algorithm then switches to another DC user and repeats the process until the target 
data rates for all DC users are satisfied. The inner loop converges because for a given DC 
user k, as L k increases, H kn for all n increases and more p kyTl become one. Hence the data 
rate R' k = J2n=i Pk,n[^og 2 (L k a kin )] + increases. On the other hand, the rates of some other 
DC users may drop due to their p fcn changing from one to zero. Nevertheless, as all the L k s 
increase, the rate of each DC user increases and hence the optimal L k 's can be approached 
iteratively. 

The complexity of the above algorithm lies in the number of iterations needed to update 
Lq in the outer loop and the number of iterations to update {L k , p k , n } for each L in the 
inner loop. Since L is searched by the bisection method, it requires log 2 (l/e) iterations 
to converge, where e is the accuracy. The empirical study in the Appendix shows that the 
exact number of iterations to find the optimal {L k , p k<n } for each L can vary for different 
values of L and different CNR realizations, but the averaged total number of iterations 
required to update the set of water levels {L , L 1; . . . , L Kl } in the whole algorithm can be 
well approximated by 0(Kf/\^N\og 2 (1/e))- Since the computational load in each iteration 
is linear in KN, the overall complexity of the proposed algorithm is 0(KjKN 1 ^ 2 log 2 (1/e)). 

D. Feasibility and Service Outage 

Similar to the delay-limited capacity problem in information theory [2], the constant-rate 
transmission of DC users considered in this paper can only be guaranteed in a probabilistic 
manner since the total transmit power is fixed and finite. The service is said to be in an 
outage if any of the basic rate requirements cannot be satisfied. Thus, the feasibility of the 
optimization problem in ©-© is directly related to the condition that P T > Pocmm, where 
-pDC,min is the minimum total power needed to support all R k s in the absence of NDC users. 
Finding Pocmm reduces to the margin adaption problem [3]. The algorithm proposed in the 
previous subsection is able to detect the service outage and obtain the outage probability 
numerically in an efficient way. Specifically, if P' T computed in Step b) of the main function 
is greater than P T when Lub is still given by the initial value set in Step a), the algorithm 



will terminate and declare an outage. 

When outage occurs, one may ignore all the K\ DC users and allocate all subcarrier 
and power resources to the K — K\ NDC users only. Alternatively, one may ignore one or 
more DC users from the user list so that at the current channel condition the resources are 
sufficient to provide the basic rates of the remaining DC users. Those ignored DC users may 
be re-scheduled for transmission at the next transmission frame by a higher layer scheduler. 
Further analysis on higher layer scheduling is out of the scope of this work. 

IV. Subcarrier and power Allocation using Dual Decomposition 

The convex relaxation technique in Section ITTT1 permits time sharing of each subcarrier. The 
system model it employs thus differs from the original OFDMA system where only mutually 
exclusive subcarrier assignment is allowed. As a result, the solution gives an upper bound on 
the achievable maximum sum-rate of all NDC users under the individual rate requirement for 
each DC user and the total transmit power constraint. Recently, it is shown in [15] that the 
duality gap of non-convex resource optimization problems in multicarrier systems is nearly 
zero if the number of subcarriers is sufficiently large. Thus, the original problem can be solved 
in the dual domain using decomposition method. Applying this result, the authors in [16] 
developed efficient algorithms to solve the weighted sum-rate maximization and weighted 
sum-power minimization problems in the downlink of multiuser OFDM systems. In this 
section we shall apply the result from [15] and solve our original problem (0 using the dual 
decomposition method. Note that the subcarrier and power allocation solution in this section 
provides a lower bound on the maximum sum-rate of all NDC users when all the target rates 
for DC users and the total power constraint are satisfied. 

Define V as a set of all non-negative -Pfc.n's for 1 < k < K and 1 < n < N such that for 
each n only one Pk, n is positive. The Lagrange dual function of the problem (0 is given by: 

g(P,[i)= max J 2 ({P kn },(3,n) (19) 
{P k ,n}ev 



max 

{Pk,n}&T> 



K N K t N 

k=Ki+l n=l k=l n=l 

K N 



\ k=l n=l 

where J 2 ({Pfc,n}, /3, Ai) is the Lagrangian and the dual variables {/3, fi} are defined in the 



same way as in ([8]). The dual optimization problem is then formulated as: 

minimize g ((3, fi) (20) 
subject to (3 >z 0,/x > 0. 

Note that the Lagrangian J 2 (. . .) is linear in /3 k and fi for fixed Pk, n , an d g((3,/j) is the 
maximum of these linear functions, so the dual problem (l20l) is convex. 

To solve the dual problem, we first decompose the dual function into iV independent 
optimization problems: 

N A'i 

g (/3, fi) = J2 9n (A n) ~ Yl ^Rk + ^Pt, (21) 

n=l k=l 

where 



g n (J3,fi)= max 



K ftTi K 

r k,n + ^2Pkr k ,n- H^2Pk,n • (22) 
_k=K 1 + l k=l k=l 

Suppose subcarrier n is assigned to user k. Using the KKT condition similar to ©, the 
optimal P£ n that maximizes the object of the max operation in (|22j) for fixed (3 and fi can 
be readily obtained and is given by dT2]) when 1 < k < K\, or (TT31) when K\ < k < K. 
Substituting (TT2l) and (TT3l) into (|22|) and comparing all the A' possible user assignments of 
this subcarrier, we obtain 



9n (/3,/x) = max 

Kfc<K 



log 2 



&k,nf3k 




/xln2 

where = for k = 1, . . . , K\ and = 1 for k — K x + 1, . . . , K. 

Once (|23l) is solved for all n's, the dual function g (/3, fi) can be obtained using (|2T|) . 
Since it is convex, a gradient-type algorithm can minimize g((3,fi) by updating {/3, //} 
simultaneously along some appropriate search directions, which is guaranteed to converge to 
the optimal solution. In general, g ((3, fi) is not differentiable, and thus its gradient does not 
exist. Nevertheless, we can resort to the subgradient derived in the following proposition. 

Proposition 1: For the dual problem (l20l) with primal defined in ©, the following is a 
subgradient of g ({3, fi) 

N 

Wk = J2 r *k,n-Rk, k = l...,K 1 , 



71=1 

N N 



av = p t -J2J2 p ^ 

k = l 71=1 

where P£ n maximizes the Lagrangian J 2 (. . .) over V at (3 and ji, and r* k n = log 2 (1 + P£ n a k , n ) 



Proof: By definition of g ((3, //) in (IT9T ): 

K N Ki / N \ / JT JV 

E r M - ^ + W ^ - EE p m 

£=1^+1 n=l fc=l \n=l / \ fe=l n=l j 

Ki / N \ / K N \ 

= ^(A/i) + E(^-^) E r M-^ p ^-EE p m • 

fc=l \n=l / \ fc=l n=l / 

Proposition 1 is hence proven by using the definition of subgradient. ■ 
With the above subgradient, both the subgradient and ellipsoid methods [17] can be used 

to update Here we choose the ellipsoid method which converges in 0((Ki + l) 2 ) 

iterations. The algorithm details can be found in [17]. The following lemma leads to a suitable 

choice of the initial ellipsoid. 

Lemma 1: The optimal dual variables {/3*,fi*} must satisfy 

0<fi*<fi max = ^- max {a k , n }, 
In 2 K^<k<K 

Kn<N 



0</^</?r x = max K, n } 

I<!<k<K 

Kn<N 



Pi 



mini^AT {«fc,n}. 

Proof: The dual variables {(3*,fi*} must satisfy the KKT conditions in order to be 
optimal. Taking the partial derivative of J2O • •) in (TT9l) with respect to P k) n results in 

In 2 (1 + a kyn P k , n ) 
if user k, for Ki < k < K, is active in subcarrier n, or 

fik&k ,n /^c\ 

i on 1 = ( 25 ) 

In 2 (1 + at k>n P k>n ) 

if user k, for 1 < k < K\, is active in subcarrier n. Since P k)U must always satisfy < 
Pk,n < Pt due to the power constraint, we obtain the upper bound /i max by letting P k , n = 
in ([24]) and the upper bound /3™ ax by substituting /i max into ([25]) and letting P fc n = P T . ■ 
Using Lemma 1, one may choose an initial ellipsoid that encloses the hyper-cuboid where 
resides, namely, E(A ,z ) = {x|(x — zo) t Aq 1 (x — z ) < 1}, where 

A = diag [ (1 + K x ) (|A max ) 2 , •••,(! + K x ) (1 + #1) (|/i max ) 2 ' 



Zf) 



1 omax 1 omax 1_ . max 

^Pl > " * ' ) 2^1 ' 2^ 



Due to duality gap, after obtaining the optimal dual variables {/3*,fi*} that minimize 
the dual function, it remains to find the optimal primal solutions {P kn } that maximize the 
Lagrangian J 2 (. . .) and satisfy all the rate and power constraints in the original problem 0. 



I / 

We can solve this by first identifying the subcarrier assignment {fi£} using (|23T) with {/3*, /i*} 
substituted and then determining the power allocation {P kn } using the results derived in 
Section III-A. 

V. Suboptimal Subcarrier and Power Allocation 

In this section we present a suboptimal allocation algorithm that has a much lower compu- 
tational cost compared with both the optimal iterative algorithm in Section III-C and the dual 
update algorithm in Section IV. The idea is to first obtain the subcarrier assignment for the 
DC users by assuming that the power is equally distributed over all subcarriers and that all the 
NDC users are absent. After that, the power distribution for each DC user over its assigned 
subcarrier set is individually refined. The purpose of the refinement is to minimize the power 
while maintaining the basic transmission rates. At last, the residual subcarriers and power are 
distributed among the NDC users using the optimal resource allocation algorithm in [4] to 
maximize the sum-rate. This algorithm is suboptimal because the subcarrier assignment for 
DC users in the first step is obtained by assuming equal power allocation. The decoupling 
of subcarrier assignment and power allocation for DC users carried out in the first two steps, 
though being suboptimal, can greatly simplify the complexity and is often used for resource 
allocation in multiuser OFDM systems such as [9], [10], [18] 

The outline of the proposed suboptimal subcarrier assignment scheme for DC users is 
presented below. 

Suboptimal Subcarrier Allocation Algorithm for DC users 

1) set R' k = 0, Q k = for all k = 1, . . . , K x and A = {1,2, ... , N} 

2) while A ^ and R' k < R k for any 1 < k < K x 

a) find k* with R' k , < R k * and R' k , - R k * < R' k - R k for all 1 < k < K\ 

b) for the found k*, find n* satisfying a k * jn » > a k * jTl for n G A 

c) for the found k* and n* , update Vt k * = Q k * U {n*}, A = A - {n*} and R' k , = 
ik+tog 2 (l + 2E f^) 

At each iteration of Step 2) in the above algorithm, the DC user whose current data rate 
R' k is the farthest away from its target rate R k will be allowed to pick a new subcarrier from 
the available subcarrier set. Preferably, the subcarrier with the highest CNR will be chosen. 

Upon acquiring £l k for 1 < k < K\, the power distribution for each DC user is then 
adjusted using the analytical solution (fl"2l) and (fl3T) . In the case where g k < \Q k \ for some k, 



the above suboptimal subcarrier assignment algorithm over-allocates subcarriers to DC user 
k. To efficiently utilize all the wireless resources, the remaining \D, k \ — g k subcarriers should 
be returned to the residual subcarrier set A, which will be distributed among the K — K\ NDC 
users. Let Pdc,t denote the actual power consumption of all the K 1 DC users. If Pdc,t is 
larger than the total power limit P T , a service outage occurs. Otherwise, the residual transmit 
power Pt — Pdc,t together with the residual subcarrier set A are allocated over the K — K\ 
NDC users. Specifically, each subcarrier in A is assigned to the NDC user with the highest 
CNR, and the power is distributed over these subcarriers in the form of water-filling (TT3T) . 
where the water level can be determined by P T — Pdc,t- 

The number of iterations involved in finding the suboptimal f^'s for k = 1, . . . , K is 
limited by N since N is the total number of subcarriers available. That is, the proposed 
suboptimal algorithm only performs a fixed number of iterations rather than iterating till 
convergence. The power allocation for given {Vt^} has explicit analytical solution as shown 
in Section III-A and its complexity is linear in KN . Therefore, the overall complexity of 
this suboptimal algorithm is only linear in K and N. 

VI. Numerical Results 

In this section we present numerical performance results of the proposed optimal and 
suboptimal resource allocation algorithms. We consider a multiuser OFDM system with N = 
64 subcarriers and K = 8 users. Therein, K\ — 4 users have DC traffic and the others 
have NDC traffic. For simplicity, we let the rate requirements of all DC users be identical 
and equal to Rqc/K\ bits/OFDM symbol, where Rbc denotes the sum of the basic rates. 
In all simulations, the channel from the base station to each user terminal is modeled by 
the HiperLan/2 channel model A [19], which is an 8-tap channel with exponential power 
delay profile, 20MHz sampling frequency and 50ns rms (root-mean- square) delay spread. 
The channels for different users are assumed to be independent. We also assume that the 
path losses from the base station to all user terminals are the same. The average channel gain 
on each subcarrier is normalized. The system total transmit SNR is defined as Pt/(N B). 
The SNR gap in the rate function (Q]) is set to 6.6 (8.2 dB) for both DC users and NDC users. 
In practice, when uncoded QAM constellation is used the SNR gap of 8.2 dB corresponds 
to a BER requirement of 10~ 5 . 

To evaluate the performance of the three proposed adaptive resource allocation algorithms, 
we also present the results for two non-adaptive schemes in comparison. In both schemes, 



the subcarrier assignment is predetermined but the power allocation for each user over its 
predetermined set of subcarriers follows the optimal approach derived in Section ITlI-AI In the 
first scheme, all the 8 users are treated equally and each is assigned 8 subcarriers. We refer 
to this scheme as Fixed Subcarrier Assignment with Optimal Power Allocation (FSA-OPA). 
In the other scheme, DC users are given higher priority than NDC users and each is assigned 
12 subcarriers, whereas each NDC user is allocated 4 subcarriers only. This scheme is called 
Fixed Subcarrier Assignment with Priority and with Optimal Power Allocation (FSAP-OPA). 
In addition, for both schemes, we let the predetermined subcarriers for each user spread over 
the entire bandwidth in a comb pattern [20]. This prevents the situation where all subcarriers 
of a user are in deep fade. 

We first compare the performance in terms of service outage behavior. Fig. [3] illustrates 
the service outage probability versus total transmit SNR when the total target transmission 
rate of DC users is Rnc = 80 bits/OFDM symbol. It is first observed that the time-sharing 
based optimal algorithm and the dual method perform almost identically. This suggests that 
two algorithms result in almost the same subcarrier assignment solution for DC users. This 
observation is expected because the duality gap vanishes when N is sufficiently large and, 
as a result, both the upper bound given by the optimal algorithm with time sharing and the 
lower bound given by the dual method approach the truly optimal solution. One can also 
see that the performance loss due to the suboptimal subcarrier assignment in the suboptimal 
algorithm is marginal. In particular, at a service outage probability of 1%, the SNR loss 
is within 0.5 dB. In addition, it is seen from Fig. [3] that the proposed adaptive algorithms 
significantly outperform the two fixed subcarrier assignment (FSA) schemes. At moderate 
and high SNR regions, the service outage probability is more than an order of magnitude 
lower. Besides, the FSA scheme with priority outperforms the one without priority as more 
subcarriers are assigned to DC users in the former. 

In Fig. HI we plot the minimum required total transmit SNR for different R DC at a given 
service outage probability of 1%. It is again observed that the optimal algorithm with time 
sharing and the dual method have almost identical performance. Therefore, only the results of 
the former will be illustrated hereafter. From Fig. [4] we observe that, for a wide range of /?dc 
that the multiuser OFDM system can support with 1% outage probability, the difference on 
the minimum required SNR between the optimal and suboptimal algorithms is consistently 
less than 0.5 dB. In particular, as -Rdc decreases, the performance of the suboptimal algorithm 
becomes closer to that of the optimal algorithm. This is expected as the suboptimality of the 
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proposed suboptimal algorithm lies only in the subcarrier allocation for DC users. If the rate 
requirement for DC users is small, the suboptimal algorithm will become nearly optimal. 
Fig. |4] also shows that, as _Rdc increases, the minimum required total SNR of the proposed 
adaptive algorithms increases at a much lower speed than that of the two FSA schemes. 

We next study the achievable transmission rates of the heterogenous multiuser OFDM 
system with the proposed adaptive subcarrier and power allocation algorithms. Fig. [5] shows 
the achievable pairs of the basic sum-rate for DC traffic -Rdc an d the average sum-rate 
for NDC traffic -Rndc at a total transmit SNR of 20 dB. The average sum-rate for NDC 
traffic -Rndc is obtained by averaging the instantaneous sum-rates of NDC users over 500 
independent channel realizations. To ensure a service outage probability of 1% or below, 
the maximum value of R DC in our simulation is set to 176 bits/OFDM symbol for the 
proposed algorithms, and to 112 and 80 bits/OFDM symbols, respectively, for FSAP-OPA 
and FSA-OPA. The maximum achievable -Rdc with an acceptable service outage probability, 
for example 1%, at a given total SNR can be obtained from Fig. |H From Fig. [5] one can 
observe that, compared with the optimal subcarrier and power allocation algorithm, the loss of 
the average achievable NDC traffic rate at a given _R DC by using the suboptimal algorithm is 
within 2% ~ 9%. On the other hand, compared with the two FSA schemes, both the proposed 
optimal and suboptimal adaptive algorithms demonstrate substantially larger achievable rate 
regions. We also observe that, at the same -Rdc, the -Rndc of FSA-OPA is larger than that 
of FSAP-OPA. This is because NDC users have fewer subcarriers in FSAP-OPA. However, 
the maximum -Rdc FSAP-OPA can support is larger than that of FSA-OPA. 

Finally, we demonstrate the multiuser diversity exploited by our algorithms. We let the 
number of DC users in the system be fixed at K\ = 4 and vary the number of NDC users 
between 4 and 16. Fig. [6] presents the average sum-rate -Rndc as a function of the number 
of NDC users at _R DC = 32 bits/OFDM symbols. Same to Fig. [51 the total transmit SNR is 
20 dB and -Rndc is obtained by averaging 500 independent channel realizations. In the two 
FSA schemes, the subcarrier allocation for DC users is the same as before, but the rest of 
the subcarriers are all allocated to one NDC user, which is selected in a round-robin fashion 
at each transmission frame. The values of -Rndc f° r the two FSA schemes remain constant 
since no multiuser diversity is achieved. On the contrary, -Rndc obtained by the proposed 
adaptive algorithms increases as the number of NDC users increases, which clearly shows 
the multiuser diversity. In particular, the achievable -Rndc of the optimal algorithm is about 
110% and 140% higher than that of the FSAP-OPA scheme when the system has 4 and 16 



NDC users, respectively. 

VII. Conclusion and Discussions 

Supporting communication services with diverse QoS requirements in future broadband 
wireless networks is crucial. This paper considered the resource allocation problem in an 
OFDM based downlink system that supports simultaneous transmission of users with DC 
traffic at constant rates and users with NDC traffic at variable rates. We investigated this 
problem from the physical layer perspective and aimed to maximize the sum-rate of NDC 
traffic while maintaining individual basic rates of DC traffic for each channel realization under 
a total power constraint. It was shown that the optimal power allocation over the subcarriers 
in such a heterogeneous system has the form of multi-level water- filling; moreover, the set of 
valid user candidates competing for each subcarrier consists of only one NDC user but all DC 
users. We converted this combinatorial problem with exponential complexity into a convex 
problem using the time-sharing technique and developed an efficient iterative algorithm with 
polynomial complexity. We also solved the original problem using dual decomposition method 
which leads to polynomial complexity as well. To further speed up the resource allocation 
and make it more suitable for practical systems, we then proposed a suboptimal algorithm 
whose computation load is only linear in the number of users and subcarriers in the system. 

The performance of our algorithms was evaluated in terms of service outage probability, 
achievable DC and NDC traffic rate pairs, and multiuser diversity. The numerical results 
showed that the convex relaxation technique with time sharing and the dual decomposition 
approach obtained almost the same solution and that the suboptimal algorithm has the 
near optimal performance. Results also demonstrated that the proposed adaptive subcarrier 
and power allocation algorithms significantly outperform the schemes with adaptive power 
allocation but fixed subcarrier assignment. 

This paper adopted the continuous rate function which greatly helped to derive the 
insights of optimal resource allocation. If discrete rates are used in practical systems, our 
algorithms can be modified accordingly. In particular, since the proposed suboptimal algo- 
rithm has near-optimal performance at significantly lower complexity, it is more desirable 
to modify the suboptimal one. For instance, one can obtain the subcarrier assignment using 
the proposed suboptimal algorithm and then apply the greedy bit loading algorithm for each 
single user as in [3]. Nevertheless, our continuous rate formulation provides the performance 
upper bound for systems with discrete rates. 



We have also assumed that the channels from the base station to all the NDC users have 
the same path loss. By symmetry, our formulation also equalizes the long-term average 
throughput among all the NDC user. To achieve fairness when their channel path losses are 
different, we can simply modify our cost function © by dividing the channel-to-noise ratio 
with the path loss. By doing so, the effective channel gains for all NDC users are normalized. 
Therefore, only the user whose current channel condition is at its peak level will be selected 
to compete with DC users for each subcarrier. This is similar to the concept of "riding on 
the channel peak" in opportunistic scheduling. 

Appendix 

Empirical Study on the convergence speed of the iterative algorithm in 

Section III-C 

Simulation settings: P T = 100, 2 NDC users, R k = 16 bits/ OFDM symbol for 1 < k < 
K x , K x 6 {1, . . . , 12}, iV e {16, 32, 64, 128}. a Kn = N\h Kn \ 2 are randomly generated with 

n modeled as complex Gaussian variables of zero mean and unit variance and independent 
for all k and n. The accuracy of bisection-searching L is set to e = 10~ 7 , and it leads to 
26 iterations in the main function throughout this simulation study. The parameter A in Step 
3b) of the inner function is set to 0.005. Note that it typically takes very few iterations in the 
while loop of Step 3c) in the inner function to find the Lk* for DC user k* that meets its rate 
requirement. Thus we choose to count the number of times it repeats for Steps 2)-3) as the 
number of iterations to update {L k , Pk, n } at a given L . Fig. [7] shows the snapshot of iterations 
to update {L k , Pk,n} at each updating step of L for three random channel realizations. The 
number of iterations varies for different L and different channel realizations and, in general, 
more iterations are needed when K x increases. To extract the rules on how the number of 
iterations change with N and K\, we plot in Fig. [8] the averaged total iterations needed to 
find the set of optimal solutions {L , Lk, Pk,n}, where each value is obtained by averaging 
over 20 independent channel realizations. For comparison we also plot the curves generated 
using the analytical expression cK\j\[N with the constant c being c = 26 x 5.1 = 132.6. It is 
observed that the analytical expression provides a very good approximation on the shape of 
the simulated curves. Therefore, we conclude that the proposed time- sharing based optimal 
subcarrier assignment algorithm converges in 0(Kf/\^N\og 2 (1/e)) iterations. 
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Fig. 1. Block diagram of a downlink multiuser OFDM system 



Z3 



II a 




Subcarriers n=l, 8 



Fig. 2. Illustration of multi-level water-filling for given subcarrier assignment in a multiuser OFDM system with 2 DC 
users, 2 NDC users and 8 subcarriers. 




Total transmit SNR [dB] 

Fig. 3. Service outage probability versus total transmit SNR at Roc = 80 bits/OFDM symbol. 




Fig. 4. Minimum required total transmit SNR versus Rue at a service outage probability of 1%. 




Fig. 5. Achievable (Rnc, Rnuc) rate pairs at a total transmit SNR of 20 dB. 
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Fig. 6. Achievable -Rndc versus the number of NDC users at a total transmit SNR of 20 dB with 4 DC users. 




Fig. 7. Iterations required to update {Lk,pk,n} at each Lo update. 




Fig. 8. Total iterations required to update {Lo, £fc, Pk.n}- 



